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ON ELECTROSTRICTION OF A GRANULAR SUPERCONDUCTOR 
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Zero-temperature field-induced polarization, supercurrent density, and the related electro striction (ES) of a granu- 
lar superconductor are calculated within a model of 3D Josephson junction arrays. Both the "bulk-modulus- driven 
ES" (the change of the sample's volume in the free energy upon the applied stress) and the "change-of-phase ES" 
(due to the stress dependence of the weak-links-induced polarization) are considered. In contrast to magnetostric- 
tion of a granular superconductor, its electroelastic behavior is predicted to be dominated by the former contribution 
for all applied fields. 
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Some attention was given recently to rather pecu- 
liar electric-field induced phenomena, either observed 
experimentally (like a substantial critical current en- 
hancement or predicted to occur (like a possibil- 
ity of magnetoelectric effect due to the Dzyaloshinski- 
Moria type coupling between an applied electric field and 
an effective magnetic field of circulating Josephson cur- 
rents |Q) in granular superconductors and attributed to 
their weak-link structure. At the same time, as compared 
to the magnetoelastic behavior of superconducting ma- 
terials (dominated either by a vortex response or 
by weak-links structure ||), their electroelastic behavior 
still remains to be properly addressed. 

In the present communication, another interesting 
pnenomenon related to the modification of the sample's 
weak-links structure in an applied electric field is dis- 
cussed. Namely, we consider a possible role of Joseph- 
son junctions in low-temperature behavior of the field- 
induced polarization and the related electroelastic prop- 
erties of granular superconductors. 

As is well-known |T(| , the change of the free energy of 
a superconductor in the presence of an external electric 
field E reads 

AF(Ei) = F(Q) - F(Ei) — V dEP(E), (1) 

Jo 

where P(E) is the electric polarization of a granular su- 
perconductor at zero temperature (see below) , V its vol- 
ume, and the internal field Ei is related to the applied 
field E via an effective dielectric constant e, namely 
Ei = E/e. When a superconductor is under the influence 
of an external (homogeneous) stress a, the above free en- 
ergy results in the associated strain component (in what 
follows, we consider only a strain component U normal 
to the applied electric field E) 

fdAF\ 



Neglecting a possible change of the effective dielectric 
constant e with the stress, Eqs.(l) and (2) give rise to the 
following two main contributions to the electrostrictive 
(ES) strains, namely 

(a) the "bulk-modulus-driven ES" due to the change 
in the free energy arising from the stress dependence of 
the sample volume 



]_dV 
V~da~ 



E, 



= T71T- dEP(E); (3) 



(b) the " change-of-phase ES" due to the stress depen- 
dence of the polarization via the Josephson junction ef- 
fective surface (see below) 



U PH = 



1 8AF\ 
V da J, 



dE 



dP(E) 
da ' 



(4) 



To proceed, we need an explicit form of the induced po- 
larization P(E). And to this end, we employ the model 
of a granular superconductor based on the well-known 
tunneling Hamiltonian (see, e.g., Ref. [|TT|| ) 



/V 



where 



with 



and 



U 



V\ da J 



(2) 



H = ^ J v I 1 ~ cos 0v (*)] ' 



^{t) = 4>i 3 (0)+uj i3 (E)t, 



-> 2e - 



4>ij(o) = 4>i- 4>j, fij = fj - Tj, 



(5) 



(6) 



(7) 



(8) 
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which describes an interaction between superconducting 
grains (with phases 0j(f)), arranged in a random three- 
dimensional (3D) lattice with coordinates fi = (xi, yi, z{). 
The grains are separated by insulating boundaries pro- 
ducing Josephson coupling with energy J^- = J. The 
system is under the influence of an external electric field 
E=(E,0,0). 

The corresponding pair polarization operator within 
the model reads [Ol 



N 



(9) 



where qi — —2erii with n% the pair number operator, and 
ri is the coordinate of the center of the grain. 

In view of Eqs.(5)-(9), and taking into account a usual 
" phase- number" commutation relation, [<^j,rtj] = iSij, 
the evolution of the polarization operator obeys the equa- 
tion of motion 



dp 1 . _ . 2e x -> 

dt = ih [p ' ] = T ^ m ^ (<)r 



(10) 



Resolving the above equation, we arrive at the following 
mean value of the field-induced polarization 

1 ■ 



P(E) = ^< p{t) > 



N 



*3 

where < ... > denotes a configurational averaging over 
the grain positions, while the bar means a temporal av- 
eraging with a characteristic time r (see below) . 

To limit ourselves with field-induced polarization ef- 
fects only, we assume that in a zero electric field, P = 0, 
and thus 0y (0) = 0. 

To obtain an explicit form of the field dependence of 
polarization, let us consider a site-type positional dis- 
order allowing for weak displacements of the grain sites 
from their positions of the original 3D lattice, i.e., within 
a radius d « ^/S (S is an effective surface of grain- 
boundary Josephson junction) the new position is cho- 
sen randomly according to the normalized (separable) 
distribution function f(r) = f (x) f (y) f (z) . It can be 
shown Jt3[ that the main qualitative results of this paper 
do not depend on the particular choice of the probabil- 
ity distribution function. Hence, assuming, for simplic- 
ity, a normalized exponential distribution law, f(x) = 
(l/d)e- x / d (where x > and / °° dxf(x) = 1), we find 

that the electric field E — (E, 0, 0) (applied along the 
x-axis) will produce a non-vanishing longitudinal (along 
x-axis) polarization vector P = (P, 0, 0) with 



P(E) = P G(E/E ), 



(12) 



, 1 / arctanz 

G{z) = - 1 

z \ z 



(13) 



Here P = 2JN/E Q V, E = h/edr, and z = E/E . 

At the same time, as is well-known ]l2| ], the super- 
current density through the Josephson junction between 
grains i and j is related to the polarization operator p as 
follows (see Eq.(10)) 



ME) 



- 1/4? 
V\dt 



(14) 



Repeating the above-discussed averaging procedure, we 
find for the change of the longitudinal part of the super- 
current density in applied electric field 



where 



j s (E) = j D(E/E ), 



D(z) = 



l + z 2 



(15) 



(16) 



where 



with j = 2eJNd/hV, and z = E/E . 

Figures 1 and 2 show the field-induced behavior of 
the normalized polarization P{E)/Pq and supercurrent 
density j s (E)/jo, calculated according to Eqs.(12) and 
(15), respectively. Notice a rather pronounced peak at 
E/E ~ 2 for both dependences. Assuming d fa lfim 
and t m 10 _16 s for an average grain size and a low- 
temperature estimate of the Josephson tunneling time 
of Cooper pairs through an insulating barrier in applied 
electric field @ in YBCO : we get E ~ 10 7 V/m for 
the estimate of the model characteristic field, which is 
very close to the typical applied field values where the 
critical current of ceramic samples was found UJ| to 
reach its maximum. Besides, assuming V ~ Nd 3 and 
taking into account that the Josephson energy in YBCO 
is J/ks — 90K , we obtain quite a reasonable estimate for 
the model characteristic critical current density jo, typ- 
ical for ceramic samples |l],|2| . Namely jo ~ 2eJ/Kd 2 ~ 
10 6 A/m 2 . Let us briefly comment on the temporal av- 
eraging (used in Eq.(ll)) and discuss the relationship 
between the characteristic time r = YijeE^d and period 
of oscillations T(E). The latter is defined via Eq.(7) as 
T{E) = 2tt/ < uo r] {E) >, where < 0Jij{E) >= 2eEd/h. 
Thus, depending on the strength of an applied electric 
field, the period of oscillations T(E) can be larger or 
smaller than the tunneling time r. In particular, high- 
field region E > Eg where the most interesting effects 
take place, is characterized by faster oscillations with the 
period T(E) < r, as compared with low-field behavior. 
In addition, we can compare r with a zero-field (and low- 
temperature) Josephson tunneling time tq ~ h/J which 
in YBCO gives tq ~ 10 _13 s. Hence, at low applied fields 
(when E <C Eq) T(E) ~ tq while in high-field regime 
(when E ^ Eq) T(E) <C t <C r . Furthermore, accord- 
ing to Eqs.(3) and (12), the explicit form of the "bulk- 
modulus-driven ES" contribution is as follows 



2 



UbMd(Ei) — UoGBMD(Ei/E ), 



(17) 
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where 

Gbmd(z) = zG(z) + In ^1 + z 2 (18) 

Here Uq = kEqPq with k — —dlnV/da being the com- 
pressibility coefficient. 

To evaluate the "change-of-phase ES" contribution 
UpH(Ei), we have to account for the stress dependence 
of the effective surface S(a) of grain-boundary Joseph- 
son junctions which was found |l5j ] experimentally to de- 
crease with the applied stress. Using the following chain 
of evident relations, 



dP 
~da~ 


_dP dG 
da da 


dP 


dPodV 
~ dV da P ° K > 


da 


dG 


dG dE 


~da~ 


dE da ' 


dE a 


dE dS 


da 


dS da' 



(19) 



and assuming that the sample volume V and the pro- 
jected area S are related in a usual way, S ~ V 2 ^ 3 , we 
obtain from Eqs.(4), (12) and (19) for the "change-of- 
phase ES" 

U PH (E i ) = ~U G PH (E i /E ), (20) 

where 

G PH (z) = zG(z) + 2G B md{z) (21) 

Figure 3 summarizes the predicted behavior of the 
two considered contributions (calculated according to 
Eqs.(17) and (20)), along with their total effect (solid 
line) on induced electrostriction. Notice that in con- 
trast to the magnetostrictive behavior of a granular su- 
perconductor (considered in Ref. fJ]), its electroelas- 
tic properties are completely dominated by the "bulk- 
modulus-driven" contribution over the whole region of 
applied fields. It would be interesting to verify the above- 
predicted behavior experimentally. 

In conclusion, a low-temperature field-induced elec- 
troelastic behavior of a granular superconductor was con- 
sidered within a model of 3D Josephson junction ar- 
rays. The "bulk-modulus-driven" contribution to the 
electrostriction (ES) of a granular superconductor (re- 
lated to the change of the sample's volume in the free 
energy upon the applied stress) was shown to dominate 
over the "change-of-phase" ES (related to the stress de- 
pendence of the weak-links-induced polarization) for all 
applied fields. 
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FIG. 1. The behavior of the induced polarization P/Po in 
applied electric field E/E a . 



FIG. 3. The "bulk-modulus-driven" Ubmd/Uq (dotted 
line), "change-of- phase" Uph/Uo (dashed line), and the total 
(Uph + Ubmd)/Uo (solid line) contributions to the induced 
electrostriction vs internal electric field Ei/Eo- 
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FIG. 2. The behavior of the induced supercurrent density 
js/jo in applied electric field E/E . 



